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ABSTRACT 

We explore a simple model for the high luminosity of SN 2006gy involving photon diffusion of 
shock-deposited thermal energy. The distinguishing property of the model is that the large "stellar" 
radius of ~160 AU required to prevent adiabatic losses is not the true stellar radius, but rather, the 
radius of an opaque, unbound circumstellar envelope, created when ~10 was ejected in the decade 
before the supernova in an eruption analogous to that of 77 Carinae. The supernova light is produced 
primarily by diffusion of thermal energy following the passage of the blast wave through this shell. 
This model differs from traditional models of supernova debris interacting with external circumstellar 
matter (CSM) in that here the shell is optically thick and the escape of radiation is delayed. We 
show that any model attempting to account for SN 2006gy's huge luminosity with radiation emitted 
by ongoing CSM interaction fails for the following basic reason: the CSM density required to achieve 
the observed luminosity makes the same circumstellar envelope opaque (r >300), forcing a thermal 
diffusion solution. In our model, the weaker CSM interaction giving rise to SN 2006gy's characteristic 
Type Iln spectrum and soft X-rays is not linked to the power source of the visual continuum; instead, 
it arises after the blast wave breaks free of the opaque shell into the surrounding wind. While a simple 
diffusion model can explain the gross properties of the early light curve of SN 2006gy, it predicts that 
the light curve must plummet rapidly at late-times, unless an additional power source is present. 

Subject headings: circumstellar matter — stars: evolution — supernovae: individual (SN2006gy) 



1. INTRODUCTION 

The recent extremely luminous supernova (SN) 2006gy 
has challenged our understanding of the deaths of mas- 
sive stars and may represent a new class of explosions. It 
was discovered by Quimby (2006) and was subsequently 
studied in detail by Ofek et al. (2007) and Smith et al. 
(2007). It is a Type Iln event, with strong narrow lines 
of H in its spectrum, although the details of its light 
curve and spectrum are unlike most other members of 
this class (e.g., Filippenko 1997). Its enormous total ra- 
diated energy of >10^^ ergs, the late (~70 days) peak 
of the light curve, and the observed constant expansion 
speed of Vexp —4000 km s~^ (Smith et al. 2007) place im- 
portant constraints on the model. SN 2006gy exploded 
in the host galaxy NGC 1260, which is a peculiar SO/Sa 
galaxy with ongoing star formation (see Smith et al. 2007 
and Ofek et al. 2007). 

The optical displays of all Type la SNe and most Type 
II SNe are believed to be dominated by the radioactive 
decay sequence ^^Ni— >^^Co— >^^Fe. If that is also the case 
for SN 2006gy, a mass M(56Ni)«10 Mq is required to 
produce the observed luminosity (Smith et al. 2007). 

Not surprisingly, SN 2006gy has prompted varying sug- 
gestions to account for its extreme energy budget. Ofek 
et al. (2007) proposed a SN la exploding in a dense H-rich 
circumstellar medium (CSM) where continual interaction 
with an external CSM was the power source for the lumi- 
nosity. They suggested that the dense CSM was probably 
the result of binary star common envelope ejection, and 
that the extreme energy budget may need to draw upon 
a super-Chandrasekhar Type la explosion or perhaps a 
massive star explosion. Smith et al. (2007) argued that 
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many different spectral properties of SN 2006gy and its 
CSM were incompatible with a Type la explosion, but 
were entirely consistent with the known properties of 
some very massive stars. To account for the luminos- 
ity with CSM interaction. Smith et al. (2007) proposed 
that the star suffered a tremendous explosive but non- 
terminal mass-loss event in the decade preceding the SN, 
analogous to the 19th century eruption of j] Carinae. 
Woosley et al. (2007) proposed a similar model, where 
the r] Car-like explosion in the decade preceding the SN 
could have been triggered by the pulsational pair insta- 
bility (different from a genuine pair instability SN) in a 
star with initial mass 110 M0. 

Interaction with an external CSM can, in principle, 
generate high luminosity as seen in other Type Iln SNe 
(e.g., Chugai et al. 2004). However, Smith et al. (2007) 
argued that this mechanism is difficult to reconcile with 
the relatively low progenitor mass-loss rate inferred from 
the weak X-rays detected by Chandra and the narrow 
Ha emission in SN 2006gy's spectrum. In other words, 
signs of CSM interaction are seen in SN 2006gy, but that 
interaction appears far too weak to power the visual con- 
tinuum. Even more troubling, the CSM interaction hy- 
pothesis cannot account for why the blast wave appar- 
ently did not decelerate while it was drained of more than 
j^qSi fjpgg (Smith et al. 2007). These inconsistencies, com- 
bined with the slow rise, low expansion speed, and high 
luminosity of SN 2006gy, led Smith et al. (2007) to ap- 
peal to the alternative of a pair-instability SN (Barkat et 
al. 1967; Bond et al. 1984; Heger & Woosley 2002) where 
the core of the star is obliterated and the light is powered 
by radioactive decay of ~10 M0 of ^^Ni. As we discuss 
below, however, it may be possible to account for the 
luminosity of SN 2006gy without placing such extreme 
demands on SN nucleosynthesis. 
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The most pressing question about SN 2006gy boils 
down to this: Do we need a pair instabihty SN to provide 
a consistent explanation for its extremely high luminosity 
and spectral properties? The pair-instability hypothesis 
requires the introduction of an exotic phenomenon; while 
the idea has been around for decades, it has never been 
observed and is only expected to have occurred in the 
early Universe (e.g., Hcger & Woosley 2002). Published 
models of the light curves from pair instability SNc do 
indeed predict high luminosities, slow rise times, long du- 
rations, and slow expansion speeds (e.g., Scannapieco et 
al. 2005) quahtativcly similar to SN 2006gy. The Hght 
curve shape depends on the assumed mass of the enve- 
lope, and most studies so far have been conducted for 
zero-metallicity stars with no mass loss. Further work 
is in progress, with various assumptions about the mass 
of the envelope (Nomoto et al. 2007; Kascn 2007; Yoimg 
2007). The progenitor of SN 2006gy is likely to have been 
a very massive star (Smith et al. 2007), but whether it 
was massive enough to suffer the pair instability remains 
an open question. 

Here we describe how a simple photon diffusion model 
can account for the high luminosity and long duration 
of SN 2006gy. Since the explosions of core collapse SNe 
typically deposit Eq ~ 10^^ ergs of energy in the SN 
envelope, an exceptional mass of ^^Ni is not necessarily 
required for SN 2006gy if an efficient mechanism could 
convert the energy of the initial SN blast into emergent 
light. In SNe from compact progenitors, the light from 
the blast wave itself normally falls short by orders of mag- 
nitude. The trapped internal radiation will cool by virtue 
of adiabatic expansion, and its energy will consequently 
be tranferred to the kinetic energy of the expanding de- 
bris. The fraction of the blast energy emerging as radi- 
ation will be roughly Ro/Rmax, where Rq is the initial 
radius of the progenitor and Rmax is the radius of the 
SN photosphere at peak light. Rq ranges from ^10~^ 
AU for the Wolf-Rayet progenitors of SNe Ib/c to ~1 
AU for the red supcrgiant progenitors of SNe II, while 
typically, Rmax ^ 100 AU. 

Falk & Arnett (1973, 1977) showed that a SN from a 
progenitor having a very extended envelope could pro- 
duce a light curve with a long timescale and high lumi- 
nosity, with very efficient conversion of blast wave energy 
into radiation. Below we consider this idea as a possible 
explanation for the light curve of SN20G6gy. While we 
argue that this mechanism may provide a plausible ex- 
planation for the early light curve, we caution that it does 
not negate the possibility that SN 2006gy is powered by 
some other mechanism at late times, such as ^^Co decay 
or continued CSM interaction. 

2. SHELL-SHOCKED MODEL 

We have in mind a model in which the SN progenitor 
(having mass Msn and radius <1 A.U.) is enshrouded by 
an opaque shell having mass Mgheii and characteristic ra- 
dius Ra ~160 A.U. The SN blast imparts Eq = lO^'^Esi 
ergs of kinetic energy. When the SN debris strike the 
shell, roughly Eq Msheii/2{MsN + M shell) of the total 
SN energy will be converted into radiation by the result- 
ing forward and reverse shocks. This radiation will be 
produced initially as X-rays, but will quickly thermalize 
to an approximate Planck spectrum with substantially 
lower temperature determined by the post-shock thermal 
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Fig. 1.— The light curve of SN2006gy from Smith et al. (2007) 
converted to Solar luminosities from M/j, assuming zero bolometric 
correction. The solid curve shows a simple light curve for a model of 
thermal diffusion from an extended stellar envelope (see text) . The 
dashed line shows the hypothetical ^®Co decay luminosity expected 
for 8 Mq of nickel. 

energy density. The radiation will escape with little adi- 
abatic loss if the characteristic diffusion timescale, Tdiff, 
is comparable to or less than the expansion timescale, 
Texp, of the ejecta. If the opacity is dominated by elec- 
tron scattering, the characteristic time for photon dif- 
fusion is Tdi//— (n ctt R?)/c^ where n is the electron 
density, (Tt=6.65x10~^^ cm^ is the Thomson scattering 
cross section, R is the initial radius, and c is the speed 
of light. For a uniform density spherical envelope with 
Solar composition, this simplifies to roughly Tdiff— 23 
days {M/R15), where M is the mass of the stellar enve- 
lope in Mq and R15 is its radius in units of 10^^ cm. The 
expansion timescale is just Texp— (AR/t;exp) for constant 
expansion speed. 

In a thermal diffusion model for SN 2006gy, we take 
the light curve peak at 70 days (Smith et al. 2007) as 
representative of both Tdiff and Texp- Prom the observed 
constant expansion speed of Vexp —4000 km (Smith 
et al. 2007), we have AR ~ VexpTexp= 2.4x10^"'' cm or 
about 160 AU. Thus, for Tdifj= T^^p, the mass of ejecta 
through which the radiation must travel is ~10 Mq. 

The visual light curve of SN 2006gy is unlike any other 
SN Iln observed to date. It is very broad, rounded, and 
smooth (Fig. 1; Smith et al. 2007). It looks Hke the hght 
curve of SN 1987A multipfied by a factor of ~250. Can 
this shape be explained by the photon diffusion model? 

In generating a simple analytic light curve, we borrow 
directly from the early study of Falk & Arnett (1973, 
1977), who studied this physical situation. In general, 
the peak due to photon diffusion becomes more lumi- 
nous, more stretched-out in time, and smoother with 
progressively larger initial radii. Following the scaling 
arguments above, we assume a constant density enve- 
lope of 10 Mq and initial radius i?=160 AU, with Tdijf= 
Texp= 70 days. 

The SN will appear suddenly at t = to, when the 
blast wave reaches the surface of the shell, delayed by 
a few weeks after the actual explosion occurs because 
during that time the blast wave was expanding through 
the opaque shell. The initial rise is basically that of a 
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constant-velocity expanding blackbody, so that L (x 
at early times (Fig. 1). It levels off near maximum as the 
photosphere recedes through the cooling ejecta, where we 
have adopted expressions from equations (40) and (44) 
of Arnett (1982). For a reasonable emitting photosphere 
temperature of re//=10'* K, the radius required at the 
light curve peak on day 70 to account for the luminos- 
ity of L=50xl0^ L0 is i?=320 AU. Working backwards 
with a constant expansion speed of 4000 km s~^ gives 
an initial radius of 160 AU at Iq. This agrees with our 
derived value for R above. 

The decline in luminosity following the peak at day 
70 is determined by thermal diffusion. The shape of the 
light curve depends on the adopted values of Tdiff and 
Texp- Arnett (1996) provides a highly simplified model in 
which the SN light is due to diffusion of internal radiation 
through a homogeneous expanding sphere. In this case, 
the light curve decline is expressed as 

L = imax e-K^=-''*+*'/2)/-rea=pTd.//] 

where t is the time in days after Iq. Adopting Tdiff = 
Texp= 70 days gives the decline rate shown in Figure 1, 
providing a satisfactory fit to light curve of SN 2006gy 
until about day 160. With a faster expansion speed and 
lower envelope mass, and hence, smaller values for both 
Texp and Tdiff, the light curve would evolve more rapidly. 
Recently, Quimby et al. (2007) advocated just such a 
model for SN 2005ap, which had an extremely high peak 
luminosity exceeding that of SN 2006gy, but faded much 
faster. 

The agreement of the early light curve of SN 2006gy 
with our simple analytic light curve in Figure 1 is surely 
an oversimplification of the actual situation. The den- 
sity profile resulting from the impact of the SN debris 
with an opaque shell will not be a homogeneous sphere, 
and the internal impact may continue over a timescale 
comparable to or greater than Tdiff- 

The total radiated energy of the thermal diffusion 
model in Figure 1 is i?rad=l-4x 10"''^ ergs, and the ki- 
netic energy for '^lO Mq expanding at 4000 km s^^ is 
Ekin=1.6xl0^^ ergs. It is mildly comforting that these 
two are nearly equal, since that is about the highest ef- 
ficiency we can expect. However, this is likely an un- 
derestimate of the kinetic energy of the SN because it is 
the mass of the accelerated pre-SN envelope (see below), 
and does not include the mass of the stellar ejecta them- 
selves, which may be comparable. Thus, the minimum 
total energy involved in SN 2006gy may be >3 x 10"''"'^ ergs, 
requiring a highly energetic explosion in any case. 

3. CAN THE RADIUS BE THAT BIG? 

The requirements of the shell-shocked model as an 
explanation for the early light curve of SN 2006gy are 
roughly the following: 1) the initial radius of the "star" 
at shock breakout must have been roughly 160 AU; 2) 
the radius of the photosphere at peak ^70 days after 
that must have been about twice as large, or ~320 AU; 
3) the mass of the extended envelope needs to have been 
about 10 M0. 

Of course, it is nonsense to invoke a progenitor star 
with an actual radius of 160 AU, since even the most 
extreme tenuously-bound red supergiants have radii 20 
times less (e.g., VY CMa; Smith et al. 2001). So the 



essential question then is how could the progenitor star 
have had such a huge initial radius? 

The most likely answer is that it was not the star's 
true hydrostatic radius, but rather, the photospheric ra- 
dius of an opaque circumstellar envelope. Smith et al. 
(2007) suggested that the progenitor of SN 2006gy may 
have suffered an explosive mass-loss event in the decade 
before the SN. This would have been analogous to events 
observed in very massive progenitor stars like rj Carinae, 
which ejected 10-15 Mq in about a decade (Smith & 
Owocki 2006; Smith et al. 2003). A shell of that mass is 
sufficient to convert 10"''^ ergs into light. High resolution 
spectra of SN 2006gy suggest that the expansion speed 
for pre-SN mass loss was of order 200 km s^^ (Smith 
et al. 2007). A mass of 10 Mq ejected at those speeds 
would require ~4xl0** ergs, a substantial energy release 
but still much less than that of a SN. 

In our shell-shocked model, the luminosity in the main 
emission peak of SN 2006gy comes from thermal diffusion 
of shock energy deposited into an optically thick ~10 Mq 
extended circumstellar envelope. Thus, the mass inferred 
from the diffusion timescale is the envelope mass, not the 
stellar mass. 

This brings up a critical distinction between our shell- 
shocked diffusion model and previous CSM-interaction 
models for Type Iln SNe. In an optically thin interac- 
tion model (see Chevalier & Fransson 1994), such as that 
invoked to explain SN 1994W by Chugai et al. (2004), a 
shock plows into a dense but optically thin CSM. Radi- 
ation emitted directly from the swept-up shell of post- 
shock gas accounts for the continuum luminosity of the 
SN, which can be considerably greater than that of nor- 
mal SNe II. For SN 2006gy, Ofek et al. (2007) suggeted 
a model where a Type la explosion expands into a dense 
H-rich CSM, and the continuous conversion of kinetic en- 
ergy into radiation powers the visual luminosity; this is 
similar to the models proposed for SN 2002ic and 2005gj 
by Hamuy et al. (2003) and Aldering et al. (2006), re- 
spectively. 

This type of CSM-interaction model cannot account 
for SN 2006gy because the extremely high CSM density 
required to power the radiated luminosity will also cause 
that same circumstellar envelope to be very opaque. It 
is clear that a circumstellar mass of order 10 Mq is re- 
quired to reprocess 10^^ ergs of kinetic energy into light 
(Smith et al. 2007; Ofek et al. 2007). If this much mass is 
distributed inside a sphere with a radius corresponding 
to the characteristic expansion radius of the SN (_R=160 
AU) then the particle density is roughly n=2xl0^^ cm^^ 
and the optical depth T=naR is greater than 300. This 
automatically forces the situation back into the thermal 
diffusion approximation described in this paper, because 
the radiation generated in the CSM interaction cannot 
escape freely. For such a large optical depth, even rela- 
tively hard (~10 keV) X-rays will be thermalized. The 
only way to avert this situation would be if the pre-shock 
CSM is fully recombined so that the opacity is lower than 
for Thomson scattering, allowing the massive envelope to 
be optically thin at visual wavelengths. Given the high 
radiative luminosity involved, this is unlikely, especially 
since high resolution spectra of the Ha line (Smith et al. 
2007) show that the pre-shock CSM is in fact ionized. 

In principle, then, the circumstellar matter ejected by 
the star before the SN event is so dense and opaque that 
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it mimics a gigantic red supcrgiant envelope with a ra- 
dius of more than 100 AU, although it is not boimd to 
the star and is probably expanding at 200 km s~^ (Smith 
et al. 2007). Thermal energy from the shock is deposited 
throughout this very large envelope, so that it can es- 
cape without suffering great adiabatic losses. This fact 
accounts for the exceptional luminosity of the SN when 
the shock breaks through the circumstellar envelope. 

Of course, after the shock escapes the pseudo- 
photosphere of the opaque shell, it will propagate into the 
lower-density wind that must surely reside outside that 
photosphere. It is this interaction that likely gives rise to 
the characteristic Type Iln spectrum of SN 2006gy, while 
the optical continuum is generated by the thermal radi- 
ation escaping from the shocked opaque envelope. This 
model may have the potential to explain many observed 
aspects of SN 2006gy that have been grounds to favor 
the pair-instability SN hypothesis over GSM interaction 
(Smith et al. 2007), such as the faint X-ray emission, 
the weak pre-shock Ha emission, the lack of any decel- 
eration in the shock, and the properties of the unusual 
broad P Cygni absorption feature in Ha. 

4. SUMMARY AND A PREDICTION 

The diffusion model for SN 2006gy that we present 
here does not require any new physical processes beyond 
those that account for core-collapse SNe. The only new 
ingredient is the hypothesis that the SN progenitor must 
have ejected an opaque shell having ~10 Mq extending 
to ~160 AU shortly before it exploded. According to 
our model, the light curve of SN 2006gy observed up 
to about day 170 comes entirely from energy deposited 
by the initial blast and does not require any source of 
radioactive energy or ongoing GSM interaction. 

However, this internal energy source cannot last. As 
the shocked shell continues to expand, the radiative diffu- 
sion time becomes shorter than the time since explosion, 
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after which the light curve will faithfully track the en- 
ergy deposition rate. Most of the kinetic energy of the 
SN must reside in debris having velocity >4000 km s~^, 
and thermal energy deposition from the impact of this 
debris with the expanding shell must terminate in <1 
yr. After that, there are only two other obvious sources 
of energy deposition to account for the light curve. 

The first is the decay of ^®Go, which would produce 
an exponential tail to the light curve having luminos- 
ity L «1.4xl0'*3 M(56Ni) exp(t/113.6d) ergs s-\ where 
M(^^Ni) is the mass of newly synthesized ^^Ni produced 
by the SN. If this is the explanation for the excess lu- 
minosity above the diffusion model after day 170, then 
about 8 M0 of "''^'Ni would be required, indicated by the 
dashed line in Figure 1. On the other hand, failure to de- 
tect a continued exponential tail in the visual light curve 
at late times does not necessarily give an upper limit 
to M(^SNi) produced in SN 2006gy. The bolometric flux 
may shift to near-infrared wavelengths as ejecta cool, and 
it is quite possible that dust will form in the shocked 
shell, reprocessing the luminosity into far infrared radi- 
ation. This cooling may also cause changes in the bolo- 
metric correction as the SN evolves with time, and may 
play a role in the tail after day ^-^160 seen in Figure 1. 

A second source could be continued shock interaction 
with the transparent GSM external to the opaque shell. 
We see clear signs of that interaction in the emission lines 
seen in the optical spectrum and in the faint X-ray emis- 
sion (Smith et al. 2007). GSM interaction could produce 
a lower-luminosity tail to the light curve of SN 2006gy 
as long as the GSM density remains high. 
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